%  TEST CASE BY HAND

it      = -1; 
keff    = 0.86784;% alt_keff;%1.039435608478710;%1.050277290361834;

%  data format
%    m1
%    g1 -->  D_1  SigR_1  nuSigF_1  xi_1  s_1->s_1 s_2->s_1 ...
%    g2 -->  D_2  SigR_2  nuSigF_2  xi_2  s_1->s_2 s_2->s_2 ...
%    ...
%    m2 ...

dat = [ 1.000   0.50   0.600  1.00  0.000  
        1.000   0.50   0.000  0.00  0.000 ];    
numg = 1;  
numm = 2;
xcm  = [   0   4 ];
xfm  = [     2   ];
ycm  = [   0   4  8];
yfm  = [     2   2 ];
mt   = [ 1 1 ];       
src  = [ 0 0];
          
BCL = 2*bc(1);
BCR = 2*bc(2);
BCB = 2*bc(3);
BCT = 2*bc(4);           

P0=[1 1]'; P1=[1 -1]';
IBSL = bc(1)*[ x1(1)*P0'-x1(2)*P1' x1(9)*P0'-x1(10)*P1' ]';
IBSR = bc(2)*[ x1(11)*P0'-x1(12)*P1' x1(3)*P0'-x1(4)*P1' ]';
IBSB = bc(3)*( x1(5)*P0 - x1(6)*P1 );
IBST = bc(4)*( x1(15)*P0 - x1(16)*P1 );

%    0.240076296817458
%    0.053576782965176
%    0.240076296817458
%   -0.053576782965176
%    0.186499513852282
%   -0.000000000000000
%    0.293653079782634
%    0.000000000000000
%    0.434801337198161
%    0.020992789338464
%    0.434801337198161
%   -0.020992789338464
%    0.413808547859697
%   -0.000000000000000
%    0.455794126536626
%    0.000000000000000
%    0.869949164114004
%    1.000000000000000

